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C\| ! Abstract 

o 

l/") ■ Statistical model of strongly anisotropic fully developed turbulence of the weakly com- 

pressible fluid is considered by means of the field theoretic renormalization group. The 
| corrections due to compressibility to the infrared form of the kinetic energy spectrum have 

been calculated in the leading order in Mach number expansion. Furthermore, in this ap- 
proximation the validity of the Kolmogorov hypothesis on the independence of dissipation 
length of velocity correlation functions in the inertial range has been proved. 
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X ■ 1 Introduction 



One of the oldest open problems in theoretical physics is that of describing fully developed tur- 
bulence on the basis of a microscopic model. The latter is usually taken to be the stochastic 
Navier-Stokes equation subject to an external random force which mimics the injection of energy 
by the large-scale modes, see, e.g. Ref. ||. The aim of the microscopic theory is to verify the 
basic principles of the celebrated Kolmogorov-Obukhov phenomenological theory |l| , study devi- 
ations from this theory and find the dependence of various Green functions (velocity correlation 
and response functions) on the times, distances, external (integral) and internal (viscous) turbu- 
lence scales. Most results are obtained within the framework of numerous semiphenomenological 
models which cannot be considered to be the basis for construction of a regular expansion in a 
certain small (at least formal) parameter, see Ref. [[J. 

One of the exceptions is provided by the renormalization group (RG) method earlier success- 
fully applied in the theory of critical behavior to explain the origin of critical scaling and calculate 
universal quantities (critical dimensions and scaling functions) in the form of e expansions, see 
Ref. HI. 
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The RG was applied to the stochastic Navier-Stokes equation in Refs. f|, [5|, ||, |7J. For ^ ne 
isotropic homogeneous turbulence of incompressible viscous fluid it allows one to prove the exis- 
tence of the infrared (IR) scale invariance with exactly known "Kolmogorov" scaling dimensions 
and the independence of the correlation functions of the viscous scale (the second Kolmogorov 
hypothesis), and calculate a number of principal constants in a reasonable agreement with the 
experiment. The detailed exposition of the RG theory of turbulence and the bibliography can be 
found in the review paper |J. 

As the model of isotropic incompressible fluid provides only a simplified description of real 
turbulent flows, it is interesting to generalize the model by taking into account anisotropy, com- 
pressibility, inhomogeneity, real geometry, and so on. In particular, in a number of papers the 
turbulence with the weak H |K], [TTJ] and strong JTJ) uniaxial anisotropy has been studied. It was 



shown that, in the three-dimensional space, the IR scaling regime characteristic of the isotropic 
case survived also if the anisotropy was included (in the language of the RG, this means that the 
corresponding fixed point remains IR stable). 



In Refs. [0, [14], [15], [UJ], the isotropic turbulence of compressible fluid was considered. The 
main difficulty is that the corresponding field theoretic model is not multiplicatively renormaliz- 
able, so that the RG technique is not directly applicable to it (for this reason, the results obtained 



in Ref. []n| cannot be considered reliable, see the discussion in Refs. |T^]) 



In Ref. ||14j| , the problem of non-renormalizability was solved in the frame of expansion 
procedure in small Mach number Ma = v c j c (where v c is the characteristic mean-square velocity 
and c is the speed of sound). In the first nontrivial order (Ma) 2 the problem was reduced to 
the calculation of scaling dimensions of certain nonlocal composite fields (composite operators 
in the language of field theory), constructed from the fields of the renormalizable model of 
incompressible fluid. 

Calculation of the scaling dimensions of composite operators is quite a cumbersome task. 
As a rule, their renormalization involves their mixing with each other, and in order to find the 
scaling dimension of a given operator, one has to consider the entire family of operators that 
admix to it in the renormalization procedure. The use of functional Schwinger equations and 
Ward identities, which express the Galilean symmetry of the model, simplifies the problem and in 



many cases allows us to find the dimensions exactly, see Refs. || [17], [18] . Using this technique 
for isotropic turbulence, the authors of |14| have calculated all the relevant scaling dimensions 
and, with the aid of these results, proved the validity of the second Kolmogorov hypothesis 
(independence of the velocity correlation function of the viscosity) in the leading order of (Ma) 2 . 
This is in agreement with the result obtained previously in |19[ within the approach based on 



the self-consistent equations. In Ref. ||15|| , this proof was generalized to all orders of the formal 
expansion in (Ma) 2 . 

It should be stressed that the stability of the Kolmogorov fixed point in the presence of 
anisotropy is obviously not a priori: the analysis of the <i-dimensional case shows that the stability 
is violated for d < 2.68 Jll], [12[ (the two-dimensional case requires special care, see Ref. p0|). The 
stability of the Kolmogorov regime is also destroyed for the anisotropic magnetohydrodynamic 
turbulence |H| and for the strongly compressible fluid [|16|] . 

In this paper, we study the effect of compressibility in the first nontrivial order of the ex- 
pansion in the small (Ma) 2 within the framework of a more realistic model of the uniaxial 
anisotropic turbulence. The anisotropy is not supposed to be small. Like in the isotropic case 



rjfl, the problem is reduced to the calculation of the scaling dimensions of a class of nonlocal 



composite operators in the model of incompressible strongly anisotropic turbulence considered in 
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However, the set of relevant operators in this case is much wider than in Ref. |T4 

EL |17L Il8| and the results obtained in 



Using the 



technique developed in || || [IT], [L8| and the results obtained in [p^| , we have found all the scal- 
ing dimensions exactly. The main result of the paper is the substantiation of the validity of the 
second Kolmogorov hypothesis mentioned above, for strongly anisotropic, weakly compressible 
developed turbulence in the first nontrivial order in the Mach number. 



2 The model 

In the stochastic theory of fully developed turbulence, the motion of a viscous fluid is described 
by the Navier-Stokes equation 

g[d t Vi + VjdjVi] = u Av { + ^didjVj - d{P + f h (1) 

the continuity equation 

d t Q + d j (gv j ) = 0, (2) 

and the equation of state V = V(g). Here d t = d/dt, di = d/dxi, Uf(x, t) are the coordinates 
of the velocity field, £>(x, t) is the density of the fluid, P(x, t) is the pressure, and v$ and v' Q are 
the molecular viscosity coefficients. Here and henceforth, summation over repeated indices is 
implied. 

Following the tradition of stochastic models of turbulence, the randomness in Eq. ([!]) is 
introduced by the large scale random force /j(x, t) with Gaussian statistics with zero mean and 
matrix of the correlation functions = (fifj), which will be specified later. 

We shall consider the weakly compressible fluid when the fields of the density and pressure 
can be written as sums of the mean values p, p and small fluctuations p, p: g = p + p, V = p + p. 
Without loss of generality, we take p = 1. Due to the smallness of the fluctuations, the equation 
of state can be taken in the adiabatic approximation: 

p = c 2 p } (3) 

where c is the adiabatic speed of sound in the turbulent medium. In the incompressibility limit 
one has c 2 = 00 or, equivalently, Ma = 0. 

Using the adiabatic relation ([|) the continuity equation (|2]) can be rewritten in the form 

1 1 

— d t p + -^diipvi) + = 0. (4) 
c c 

For c 2 = 00 the density becomes a constant, the velocity becomes transversal (diVi = 0), and we 
return to the case of incompressible fluid. 

The velocity field Vi can be expressed in the form Vi = + t>J, where vj- = P^ v j is the 

transversal part satisfying the condition divj- = 0, and vj = PfjVj is the longitudinal part. The 
longitudinal P" and the transversal P 1 - projectors in the wave-vector (k) space have the form 
P}j = kikj/k 2 and P^j = Sik — P}j, respectively (k =| k |). 

Since the velocity field has to become transversal in the incompressibility limit, its longitudinal 
part v" has to be proportional to the inverse square of the sound speed, |v|" ~ c -2 . Then from 
the Navier-Stokes equation it follows that |f |" ~ c~ 2 for the longitudinal part of the random 
force. Hence, c~ 2 can be treated as a small formal parameter, and the compressibility corrections 
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to the transversal part of the velocity field can be studied within the expansion in c 2 (or in 
Ma 2 ). 

In the first order in c~ 2 the continuity equation takes the form 

hdt+vid^p + div} = 0. (5) 

In the leading approximation in c~ 2 (corresponding to the incompressible fluid) the Navier-Stokes 
equation ([!]) gives the well-known relation between the pressure p and the transversal velocity 

Ap = -didjV^vf. (6) 

The last two equations allow us to express the pressure and the longitudinal part of the 
velocity via the transversal part v 

v\ = -^A-'Vtp, p = -A^didjvj-vj; (7) 

where Vt = d t + v^di is the Lagrangian derivative for the transversal part of the velocity and 
A -1 is the Green function for the Laplace operator. In the field theory the quantities like the 
right-hand sides of Eqs. (0) are termed "composite operators" . 

Operating with the transversal projector P onto Eq. (H) and using relations (0), we arrive 
at the closed equation for the transversal part of the velocity (and, therefore, for all its statistical 
moments) in which the compressibility is taken into account to the order of c -2 , or, equivalently, 
Ma 2 

d t vi = v Avi - Pjjtfd.vf] - Pjjtfd.v) + v\d s vf] - c- 2 u,P^[pvf} + ft- (8) 

To simplify the notation, we shall write Vi instead of vf- in what follows. 

The positively definite (dx d)— square matrix of the pair correlation functions of the random 
force ft will be taken in the form (see, e.g., 0, §]) 

(//(x, t) ft(0, 0)) = e D js (x, t) = 5(t) s J^- d D? s (k) exp [ikx] , (9) 

A,(k) = fc 4 - rf -^(k) (10) 

(we recall that (/• ) = 0). We see that temporal correlations of ff have the character of white 
noise, while the spatial falloff of the correlations is controlled by the parameter e and space 
dimension d. The functions (^) are translation invariant and for e = 2 become scale invariant 
when the amplitude Eq acquires the dimension of the energy dissipation rate, e, see Ref. 0. The 
value e = 2 is physically most acceptable, since it represents the assumption that random force 
acts at very large scales, which substitutes for effect of boundary conditions. For simplicity, we 
use the force correlation function (|^) without the usual infrared regularization. In this case, for 
e = 2 the function (|10|) with the proper choice of the amplitude Eq in Eq. ® can be considered as 
a power-like model of the "ideal" pumping function 5(k), see ||. The justification of this choice 
as well as the discussion of the central problem of the e— expansion, i.e. the continuation from 
to e = 2 have been thoroughly discussed in Ref. p4 



The ratio So/v 3 = g plays the role of a bare coupling constant, i.e. the expansion parameter 
in the nonlinearity (v<9)v in the non-renormalized perturbation theory. In the limit e — > the 
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constant go becomes dimensionless, the diagrams of the Green functions become divergent in the 
ultraviolet (UV) region of the wave-vector space, and the problem of eliminating these diver- 
gences emerges. In the field theory this problem is solved by the well-known UV renormalization 



procedure, see, e.g., 



In this paper we consider the uniaxial anisotropic turbulence. The transverse projector V 
for the correlation matrix (|T0|) is defined by the relations || 11], |12jl : 



Vj- S (k) = (1 + c^P^lc) + a 2 i^(k) , (11) 
j£(k) = Sj.-P^k), Pj s (k) = kj k s k' 2 , 

i? is (k) = (n,- - K - Cfe^X) , tk = (kn)k-\ (12) 

where the unit vector n yields the direction of the anisotropy axis and ax, «2 are free amplitudes. 
These amplitudes are not considered small in the present analysis, however, restrictions to their 
values «i > —1, «2 > follow from positive definitness of the matrix (|9|). For nozero olx,oli the 
random forcing describes differences in energy injection in the prefferred direction and directions 
perpendicular to it with the subsequent generation of anisotropic structures in large-scale eddies. 



3 Field theoretic formulation and the RG equation 



As in the critical dynamics f2~2| , the stochastic problem (|8D,(P),(|10|) is mapped to a quantum- field 



model, which is determined by an effective De Dominicis-Janssen "action" S(v, v') constructed 
on the basis of the original stochastic model. This action is a functional of the transversal velocity 
v and an independent transverse auxiliary field v'. 

In this approach, the generating functional G of the velocity correlation and response functions 
is the functional integral 

G(A, A') = j Dv DV det M(v) exp [S(v, v') + Av + A'v'] , (13) 

with the effective action 

5(v, v') = -go^v'Dv' + v'[-d t v + z/ Av - (vd)v - (v"d)v - (vd)v 11 - C - 2 i/ Avp], (14) 

where A, A' are the source fields, which are equivalent to regular external forces. Here, the 
required integrations over the spacetime arguments of the fields and sums over discrete indices 
are implicitly assumed. 

The Jacobian det M in Eq. ( |13|) ensures the cancellation of all the diagrams containing the 
self-contracted bare propagator (vv'), which arise along with other diagrams from the rules of 
the Feynman diagrammatic technique for the action (fUD, but do not arise in the construction of 
diagrams by direct iteration of the stochastic equation (§). Following |22|, |5], [|, we simply define 
these superfluos diagrams as zero, and simultaneously set det M — 1 in Eq. (p~3"D . We note that in 



our model such a definition is nontrivial because the interaction in (14) involves the derivatives 
with respect to the time variable. Nevertheless, this definition is feasible, as it has been shown 
in Ref. [14[ using isotropic turbulence as an example. As a result, we arrive at a standard field 



theoretic model with action (13), and the standard renormalization theory is applicable to it. 
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The action fll4|) is not renormalized and the corresponding Green functions of the fields v, v' 
contain UV divergences for e — > 0. In order to analyze them, we rewrite the action ( |14|) as a sum 
S = S 1 + S c : 

S 7 (v,v') = ^ 3 vW + v'[-9 t v + . Av-Hv], (15) 



S 



c, 



v,v 



(16) 



where a i = — c 2 ,ao2 = c 2 z/ , and the composite operators Fi,F 2 , according to Eq. (Uh and 
using the relation d^j = djvf, can be represented in the form 



In the limit c~ 2 



,11 



v'^djVi - diVj)diA 1 V t A 1 d i d j v i v j , F 2 = v[(Avi)A 1 d i d j v i v j 



(17) 



0, the action fll5|) describes the incompressible anisotropic turbulence. Renor- 

It was shown that in order to ensure 



malization of this model has been considered in 11, LL 



the multiplicative renormalizability, the model has to be extended by adding certain anisotropic 
dissipative terms with new viscosity coefficients i^oXoi, i — 1,2,3, where the dimensionless pa- 
rameters Xoi describe the relative impact of the different anisotropic structures on the viscous 
dissipation and play the role of additional coupling constants. 

The renormalized action corresponding to the original non-renormalized functional (O) is of 
the form 



S 



= -gu 3 fi 2t VDV + v'[-<9 t v + vZ v A-w + z/Z,xiZ Xl nA(vn) + 
+ uZ vX 2Z X2 {nd)\ + vZ uX 3Z x Mnd) 2 (vn) - (vd)v\. 



(18) 



Here, the renormalization mass fi is an additional arbitrary parameter of the renormalized theory, 
the renormalized parameters g, u, Xi are related to their bare (unrenormalized) counterparts by 



the multiplicative renormalization formulae, |TJ 

go = gZ g fi 2e , uq = uZ u , xoi = XiZ- 



'-3 



Xr 



(19) 



The renormalization constants Z are calculated within the perturbation theory. In the minimal 
subtraction scheme they have the form U Z— 1+ only poles in e" and cancel all the UV divergences 
in the correlation functions of the primary fields in the model (0). The last relation in ([19|) follows 
from the absence of the constants Z in the first and last terms of the action (ffH). 

To determine the dependence of the renormalized correlation functions on the parameters a i 
and a 02 after the term ( |T6"D has been added to the action, let us consider the pair correlation 
function for the incompressible isotropic case (the detailed discussion can be found, e.g., in Refs. 

<! k Gf m (k)exp[ik( Xl -x 2 )]. (20) 



(21) 



The RG equation for the trace of its space Fourier transform G R (k) = G^Qs.) is 















G R (k) = . 



where the j3 g function and the anomalous dimension 7^ are expressed via the renormalization 
constant Z v 

P a = -g(2e - 3 7 ,), 7, = V^lnZv . (22) 
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Here £> M denotes the operation fid/dfi taken at fixed values of all the bare parameters. 
The solution of the RG equations along with the dimensionality considerations gives 



G R (k) = u 2 (s)k 2 - d R(g(s)) 



2-d 



k 

S = -, 



(23) 



where R is a "scaling function" of the invariant charge g(s), the effective variable satisfying the 
equations 

sf s =f3 g (g(s)), g\ s=1 = g. (24) 
The second effective variable, the invariant viscosity v(s), satisfies the equations 



s^T = ~lv (g(s)) ,u\ s=1 
as 



(25) 



From the solution of equations (2H|) it follows that g(s) — > g* for s — > 0, where g* is an infrared 
stable fixed point of the RG equations, i.e. the root of the equation f3 g = with the positive 
value of the correction exponent uj = dj3 g /dg. 
The solution of equation is 



u(s) = z/exp 



f3 g (x) 



(26) 



From (|26| ) along with Eqs. 



(HD it follows that 

1/3 



VK8] 
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gs 2 \ 



1/3 



(27) 



For the spectrum of kinetic energy E(k) ~ k d ~ 1 G R (k) in the asymptotic region s — ► we 
obtain from Eqs. (|23| ) and ([27]) the expression ^(/c) ~ E 2 / 3 k~ 5 / 3 , which is independent of the 
viscosity u and corresponds to the Kolmogorov value of the exponent. 

When the anisotropic case is studied [action (0)], the new terms j3 x .dG R (k)/dxj are ap- 
pended to the RG equation (21) . The new (3 functions and the anomalous dimensions 7 Xi corre- 
sponding to the new dimensionless parameters Xi 

Pxi = -XiJxn Ixi = V^\nZ Xi (28) 

are expressed via the renormalization constants Z Xi in the action The additional invariant 
variables x( s ) satisfy equations like Eq. (0). In Ref. [I2| it has been shown that those equations 
have an IR stable fixed point g(s),x( s ) ^ 9*iX*i m which all the eigenvalues of the matrix of 
the correction exponents 



dp, 



a, 



ij 



9i=g * , 9i = 9,XuX2,X3, hJ 



0,1,2,3 



(29) 



(to be precise, their real parts) are positive, i.e. the Kolmogorov asymptotic regime conserves 
the stability against the strong anisotropy. 

The problem becomes more involved if the compressibility is taken into account. Let us sup- 
pose that we have managed to renormalize the action (fL6[). Then, the new terms 7 ai P ai G /? (k) 
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appear in the RG equation ([H]), where 7 ai are the anomalous dimensions of the renormalized 
parameters Oj. In contrast with the parameters \i the renormalized counterparts of the parame- 
ters aoi and ao2 have nonzero dimensions and, therefore, the scaling function R depends on the 
effective dimensionless variables 

TTj" = k 2 T7 2 a 1 , = A; 2 z/a 2 • (30) 

The effective variables ai(s) and d 2 (s) satisfy equations like Eqs. (^). In the infrared asymp- 
totic region k — > they take on the form d{ ~ fc _7 % and the infrared asymptotic form of the 
dimensionless arguments (|30D is given by the expressions Ui ~ k~~ Aa i with the scaling dimensions 
A Qi (for more details see, e.g., |25L |2B|). In the linear approximation with respect to the small 



parameters aoianda^ in the functional ([13]) the leading correction to the scaling function R takes 
on the form (1 + const ■ k~ max ) where A max is the maximal dimension among A a .. 

Therefore, the investigation of the dependence of the kinetic energy spectrum on the com- 
pressibility is related to the calculation of the scaling dimensions A a% which, as we shall shown 
below, can be expressed via scaling dimensions of the composite operators F 1]2 entering into the 
action (1161) - 



4 Renormalization and scaling dimensions of the compos- 
ite operators 

The addition of the term (|i~6|) involving the operators F\ , F 2 (|17D to the action (pl|) gives rise 
to new UV divergences (poles in e) in the correlation functions. According to the generic rules, 
all the composite operators with the same canonical (naive) dimensions and tensor structure 
can be mixed in the renormalization procedure, i.e. an UV finite renormalized operator F R 
has the form F R = F + counterterms , where the contribution of the counter terms is a linear 
combination of F itself and other unrenormalized operators that "admix" to F. Therefore, to 
perform renormalization of the operators FiandF 2 , one has to consider a wider family of operators 
Fi which admix to F\, F 2 . 

The renormalized operators Ff 1 are related to their non-renormalized counterparts Fi by the 
well-known matrix formulae of multiplicative renormalization, see, e.g. Refs. || || 

Fi = ZijFj 1 , (31) 

where is the matrix of the renormalization constants. In the minimal subtraction scheme 
its diagonal elements have the form 1+ poles in e while the non-diagonal elements contain only 
poles. From the matrix one calculates the matrix of anomalous dimensions 7^ = Z^'D^Z^ 
and the matrix of scaling dimensions for the set of operators 

A{ = £§ + 70. (32) 

The contribution Df- = [d F — -f u dp]ij is expressed via the anomalous dimension of the viscosity 
(p^), and the total dp and frequency d F canonical dimensions of the operator F || |8|], which are 
equal to the sums of corresponding dimensions of the fields and derivatives that constitute F. 

The total canonical dimensions of the fields and parameters of the model are found from the 
requirement that all the terms of the action ([14]) be dimensionless, see || §: d t = d w = —2, 



d v = 1, d v > = d — 1 (d x = dk = — 1 by definition). From these dimensions we then obtain the 
canonical dimensions of the operators F±, F 2 equal to dp 1 = dF 2 = d + 4. We also note, that these 
operators are Galilean invariant, scalar and nonlocal. 

Let F = {Fi} be a system of composite operators closed with respect to renormalization. 
The equation e^F/* = a 0i F (the summation over the subscript i is implied) can be regarded as 
a definition of the renormalized sources a = {a^}, which for the usual renormalization formulae 
do = aZ a , F = ZpF R leads to the relations Z a = Zp 1 for the renormalization constants and 
If = —la for the corresponding anomalous dimensions. The requirement that the terms 



J dxaF R = J dxa^F , x = x, t 



be dimensionless then gives the "shadow relations" for the canonical and scaling dimensions of 
the operators Fi and sources a« 

d k a + d F = d, d£ + d£ = l, A a + A F — d + A w . (33) 



Due to Eqs. (|3"3"D, the problem of finding the maximal dimension A ai for the sources corresponding 
to the operators F\ i2 i n the action (|16|) is equivalent to the calculation of the minimal scaling 
dimension A^ associated with the operators F\ i2 and all the operators that admix to them in 
renormalization. 



According to the general theory of renormalization, see, e.g. fl21| , counter terms in a field 
theory with a local interaction are also local. Therefore, the renormalization of the nonlocal 
operators Fi, F 2 is reduced to that of their local blocks (see below) and to the admixture of the 
local operators (i.e. monomials constructed of the fields and their derivatives at the same point 
x, t) with the same canonical dimension and symmetry (Galilean invariant scalars). These local 
operators in our case are the following: F = dv'dvdv, dv'Vtvdv, dv'd 3 v, n 2 dv'dvdv, n 2 dv'Vtdv, 
n 2 dv'd 3 v, n dv'dvdv , n dv'Vtdv, n dv'd 3 v, n 6 dv'dvdv, n 6 dv'Vtdv, n 6 dv'd 3 v. The notation is 
symbolic and it implies all possible contractions of the vector indices of the fields v', v, derivative 
d and unit vector n. This set of operators is closed with respect to renormalization because the 
nonlocal operators Fi and F 2 cannot admix to them. The first three types of the operators F 



have been considered in [ FJ] . It was shown that they did not affect the scaling dimensions of the 
nonlocal operators Fi, F 2 due to the fact that the corresponding renormalization matrix was 
block-triangular. This feature of the renormalization matrix persists also in the other operators 
F, which contain the vector n, so that they also do not affect the scaling dimensions of Fi, F 2 . 
In contrast with the local operators F, they contain additional factors of A _1 dv which have zero 
canonical dimension and negative scaling dimension —4/3 at e = 2 (we recall that the scaling 
dimension of the field v equals to —1/3 at e = 2, see || |J). Therefore, the scaling dimensions of 
the operators F are greater than the dimensions of the nonlocal operators F\, F 2 , and the leading 
contribution to the IR asymptotic form of the spectrum is determined by the contributions of 
Fi and F 2 . We note that due to renormalization, scaling dimension of an operator F does not 
coincide in general with a naive sum of scaling dimensions of the fields and derivatives entering 
into F. But, for the incompressible case, the hypothesis that scaling dimension of a nonlocal 
operator is the sum of scaling dimensions of its local parts and of the factors of type A _1 <9t> has 
been confirmed in |25| by the explicit one-loop calculation of the scaling dimensions related to 
the local operators with the canonical dimension d + 4, and we also accept it in what follows. 

As result, we obtain that the scaling dimensions of F\ and F 2 are determined by their own 
renormalization. The latter is reduced to the renormalization of the local blocks entering into Fi 
and F 2 . 
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Let us denote the field v by the solid line, v' by the oriented solid line, and the operator 
A -1 by the wave line. The derivative with respect to coordinate is denoted by a slash, and 
the derivative with respect to time by a cross. Graphical representation of the operators ([17]) 
is depicted in Fig. 1, where the vector indices are omitted and the operator, containing the full 
time derivative Vf, is represented as a sum of the first two diagrams. 

The contribution of the last operator from Fig. 1 to the correlation function (v'vvv) is 
depicted in Fig. 2. The shadowed rectangle denotes an arbitrary one-particle irreducible diagram 
with fixed external legs. One can show that the triangular subdiagram contains UV divergence 
and its elimination requires the renormalization of the local block of the nonlocal operator under 
consideration. Thus, for the complete renormalization of the operators fll?!) it is sufficient to 
study the renormalization of all their local blocks. 

The operator Fx consists of two nonlocal factors A" 1 , the full derivative Vt, and two local 
blocks 

G i = v 'i(djVi - diVj), G 2 = didjViVj, (34) 
while F 2 contains one factor A -1 , the operator G 2 , and the local block 

G 3 = vlAvi. (35) 

The scaling dimensions of the operators ( |i~7D are equal to the sums of the scaling dimensions 
of the above factors, among them only the dimensions of the local blocks ( p4|) and ( p5|) require 
nontrivial calculation. In order to find them one has to study the renormalization of the complete 
set of the operators that admix to Gi in renormalization. This set is rather big because of the 
anisotropy and the canonical dimension of Gi is high (dp = 7 for d — 3). To simplify the analysis, 
we shall use some general rules for the operator mixing. Their proof and other examples can be 



found, e.g., in Refs. H, R [T7|, IT 



(a) In the action ([18|) the derivative in the interaction term can be moved onto the auxiliary field 

v' using integration by parts: v^VjdjVi = —(djv'^ViVj. Therefore, the derivative d appears 
as an external factor for each external leg of the field v' for any one-particle irreducible 
diagram, and the corresponding counter term contains the factor &v' . 

(b) Only Galilean invariant operators can admix to an invariant operator in the renormalization 

procedure. 

(c) Let some operator G has the form of a total derivative of some other operator [G], G = d[G]. 

In this case, the scaling dimension of G is simply given by the relation Aq = 1 + Ar^i. 

(d) All the one-particle irreducible diagrams, containing closed circuits of the retarded propa- 

gators (vV), vanish. 

We denote by G or [G] the full sets of operators that can mix with a given G or [G] in 
renor malizat ion . 

According to the item (c), instead of the operator G 2 from (p5Sf) it is sufficient to study the 
renormalization of the operator [G 2 ] = v(Vj. Due to the transversality of the field the only 
operators that can admix to [G 2 ] have the form [G 2 ] = nkfiidiVk5ij, nknidiVkiiirij. Their scaling 
dimensions are equal to Ar^i = 1 + A v . The scaling dimensions of the fields v, v' and the time 
have the form (see, e.g. ||) 

A„ = 1 - 2e/3, A v , = d - 1 + 2e/3, A t = -2 + 2e/3. (36) 
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We then obtain Ar G i = 2 — 2e/3, which gives Ar^i = 2/3 at e = 2. Since the operator [G 2 ] itself 
is not renormalized, for the scaling dimension of G 2 we obtain (item (c)): A G2 = 2 + Ar G2 j= 
2 + 2A„ = 4 - 4e/3, which gives A Ga = 4/3 at e = 2. 

The operator Gi consists of two terms: G\ = Gn — Gi 2 . The term G 12 is rewritten in the 
form G12 = diiy'jVj)] it is then sufficient to consider the operator [G 12 ] = v[vj (item (c)). It can 
mix with the following operators: d(v'p riiUidiVp 5ijn k nidiv' k , and riinjn k nidiv' k . They all are UV 
finite and their critical dimensions are simply given by 1 + A„/ = d + 2e/3, i.e., 13/3 at d — 3 
and e = 2. The diagonal element of the matrix Z^ of the above set of operators equals 1 (item 
(a)) and, as in the case of the set associated with G 2 , this matrix is triangular. It then follows 
that A Gl2 = 1 + A[ Gl2 ] = 1 + A v + Ay/ = d + 1, which gives A Gl2 = 4 at d — 3. 

The operator Gn does not admix to itself due to the item (a). Owing to the Galilean 
invariance, it does not mix with the operators of the same tensor structure which involve the 
vector n (item (b)). Furthermore, it does not mix with the invariant operators rijVtn k v k (item 
(a)) and riiVtdjVi (item (d)). The set Gn, which can admix to Gn, includes the operators Av'j, 
rijAv' s n s , n s d s nidiVj, djn s d a v[ni, and n^n k d k n s d s v[rii. All these operators are UV finite and their 
critical dimensions are equal to 2 + A v > = 16/3. Like the case of the operators G 2 and Gi 2 , these 
operators do not affect the critical dimension of Gn = v[djVi. Since the latter is UV finite, its 
critical dimension is given by A Gll = 1 + A v > + A^ = 4. 

Now let us turn to the last operator G3 from fl3~5|). The invariant operators t^Vtfj and 
v'jHiV 'tVjUj do not admix to G3 due to item (a). Therefore, we are left with the three types of 
operators 

{G 31 } = H(nd)\, (nv')A(nv), (nv')(nd) 2 (nv)}, 
{G 32 } = {dtiv&vs), dttv'dvi), (n<9)K(n<9H, (nd)[{Vd)(nv)\, 
^[(nv')^(nv)], «9 z [(nv')(na>], (n0)[(nv')(nd)(nv)]}, 
{G 33 } = {(nd)A(nv'), (nd)(n«9) 2 (nv')}. (37) 

The operators {G33} do not affect the scaling dimensions of {G31} and {G3 2 } (item (c)), they 
are UV finite and their dimensions are equal to A| Ga3 | = 19/3 at d — 3 and e = 2. The operators 

{G 32 } do not affect {G31} (item (c)), they are also finite (like Gi 2 ), and their scaling dimensions 
are equal to A{G 32 } = 5. 

Thus, we need to renormalize the remaining set that includes the operators G3 and G31. They 
are renormalized with mixing, and the corresponding matrix Zij is nontrivial. In isotropic case 
the renormalization constant of G3 is expressed via the known renormalization constant Z u in 
the action ([18]) and, therefore, the scaling dimension Aq 3 is related to the known function 7^ 



14j| . In the presence of anisotropy the situation becomes more complicated. However, even in 
this case it turns out possible to express the matrix Z^ in terms of the known renormalization 
constants Z v and Z Xi from the action ([18]). 

Consider the generating functional ([T3|) with det M = 1 and the renormalized action (|T8"P. It is 
UV finite and, therefore, its derivative with respect to the renormalized parameters e = {g, Xi, v} 
(they are the generating functionals of the composite operators d e S) are also UV finite, as well 
as the operators d e S themselves. 

The functional G(A, A') satisfies the RG equation 

D* G G(A, A') = 0, V RG = fr± - 7 „„|- + (3 g ^- + & A] (38) 
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with the functions (3 g and f3 Xi defined in ( p2]) and (^). Let us define the matrix u>ik by the 
relation 

w« = -ft^, (39) 

where $ = {(?, Xi} ( we recall that 7 5 = —37,,). Using ( p2|) and we readily find that at the 
fixed point g* the matrix (j39|) coincides with the matrix of correction exponents defined in 



We define the commutator of two differential operators T>i, T>j in a standard way, [Di,T>j] = 
ViVj — VjVi. The commutators of the operators T>ng, T> v = = ud/du and <9 9i = d/dgi are 
of the form 

Pro, V v \ = 0, [Phg, c> gi ] = ^[(Sio^-P, - flJ ] - fo^,. (40) 

Differentiation of the RG equation ([38|) with respect to v and gi along with the commutation 
relations fl40| ) gives 

V RG d 9i G = oo tJ [5 t0 ^V u - d 9j ]G - 5 l0 ^d g G, V RG V U G = 0. (41) 
&9 9 

The fact that the operators T> RG and T> v are commutative allows the left-hand side of the first 
equation in (|4l"l) to be rewritten in the form 



V RG d 9i G = V RG [d gi - 5 iQ ±-V v \G - S m ^V v G. 



Using this relation, equation fl4l|) is rewritten as 

V RG X { = -UijXj - 5 i0 ^-X , i,j = 0, 1, 2, 3, 4 (42) 
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which, at the fixed point g* ^ along with f3 g = 0, gives 

V RG Xi = -UijXj . (43) 
This is nothing else than the scaling equation for the quantities 

X i = (d gi -(3g)- 1 8 i0 V v )G(A,A'), 
and Uij is the matrix of their anomalous dimensions. Its eigenvalues are positive (it follows 



from the IR stability of the fixed point, see |T2[)- According to (p9|) , it is expressed via the 
renormalization constants Z of the action ([18]) calculated in the one-loop approximation in Ref. 



0- 

Using the explicit form of the generating functional (|T3"D, the quantities Xi are explicitly 
expressed via the derivatives of the renormalized action ( |18D with respect to the parameters g, 
v, and Xi 

X i = J DvDv'Xi exp[S(v, v') + Av + AV] = 
= J DvDv'[d gi S - {?>gy l 5 i0 V v S} exp[S(v, v') + Av + AV]. (44) 
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Therefore, the quantities Xi represent the generating functionals of the correlation functions that 
involve not only primary fields v, v' but also renormalized composite operators Xj. 
Performing the differentiation in (fHj) explicitly one obtains 

(45) 



Xi = e i0 v'Av + ea(v'n)A(vn) + ej 2 v'(n<9) 2 vei 3 (v'n)(n<9) 2 (vn 
where the coefficients e are expressed via the renormalization constants from ( 



e o 


= v [d g Z v 


z v \ 
3g) 


\g=g* 


> e 0i - 


= "Xi 


d 9 


(z v z x . 


eio 


= fd Xl Z u 


9=9* ' 


eii 


= "Xid Xl 


(Z u Z Xi ) 


9=9* ' 


e 2 o 


= vd X2 Z v 


9=9* ' 




= "Xid X2 


(Z U Z X .) 


9=9* ' 


e30 


= vd X3 Z u 


9=9* i 


Z3i 


= "Xid X3 


(z„z Xi ) 


9=9* ■ 



ZyZ Xi 

3^7 



\9=9 i 



i = 1,2,3 



(46) 

It is obvious from Eqs. fl46l ) that are given by linear combinations of the operators G3 
and {G31}, and the matrix u (^) determines their anomalous dimensions. The eigenvalues u>i 
of the matrix uj have been calculated in Ref. |l!| in the first order of the e expansion. All the 
real parts of these eigenvalues are positive (two of the eigenvalues are complex). We calculate 
from Eq. (|32]) the scaling dimensions of the operators G3 and {G31} Ag 3 = 13/3 + uj (uj = uq), 
A 



rg 3l j. = 13/3 + uji for i = 1,2,3. From the results of Ref. 11211 it follows that the exponent 
uj is smaller than each of the eigenvalues uji, so that the main contribution of the operators in 
consideration to the IR asymptotic behaviour of the kinetic energy spectrum is given by the 
operator G 3 . 

Finally, from Eqs. (|T7|) , fl34|), (|35|) and "shadow relation" ([33|) we obtain the scaling dimensions 
for the original composite operators F\, F 2 and the corresponding parameters a\anda 2 



A 



ax 



d + 4 - 2e, 
4e/3 - 2, 



A 



t2 



d + A-2e + uj, 



A 



"12 



4e/3 



UJ. 



(47) 
(48) 



For d = 3 and e = 2 this gives A ai = 2/3, A a2 = 2/3 — uj (A a2 = —10/3 in the first order in e) 
Since the parameter ax is not renormalized (see above), we have ai = a\ = aoi = — c _ 
which along with Eqs. (|30| ) and (B^) in the IR asymptotic region for the effective variable u\ 



gives: Ui 



,-2-1 



£ i/3^,-2/3_ Tj s i n g well-known relation £ ~ v^/L (where e is the mean 
dissipation rate, v c is the mean-square velocity field, and L is the outer scale of turbulence) the 
expression for u\ can be rewritten as 



U-, 



(Ma) 2 (fcL)- 2/3 . (49) 

In a similar way, one can find the k dependence of the variable u\ at e = 2. From the relation 
u 2 ~ u\(klY (see [3(]), where Z = £ _1// Vq' /4 is the Kolmogorov dissipation length, one obtains 

u* 2 ~ (Ma) 2 (E)- 2/3 (fc!) w , 



(50) 



and >> «2 (^ > 0) in the inertial range kl << 1. Therefore, the leading contribution to the 
small k behavior of the scaling function R from Eq.(|2"3"|) is given by the term with the variable 



In the linear approximation in the Mach number, the leading correction to the kinetic energy 



spectrum is of the form 



E(k) ~ e 2/3 k~ 5/3 \l + A(Ma) 2 (kL) 



-2/3 



(51) 
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where A is a numerical factor. This correction is independent of the viscosity coefficient vq which 
proves the validity of the second Kolmogorov hypothesis. The contribution of u* 2 gives rise to 
a uq dependent term, but in the inertial range it only determines a vanishing correction. For 
Ma << 1 the correction is rather small because in the inertial range one has (kL)- 2 ' ?J < 1. In 
contrast with the isotropic model, the amplitude factor in ( [51] ) and the coefficient A depend on 
the anisotropy parameters. 



5 Conclusion 

We have shown that in the statistical model of the fully developed turbulence in the presence 
of uniaxial anisotropy the kinetic energy spectrum in the inertial range is independent of the 
viscosity coefficient (i.e. the second Kolmogorov hypothesis holds) in the leading approximation 
in the Mach number. 

In this paper, we have dealt only with the dependence on the UV scale (or, equivalently, on 
the viscosity coefficient) and have not discussed the dependence on the integral scale L. The 
RG approach along with the operator product expansion is also applicable to this problem. The 
most singular L dependence is revealed by the different-time velocity correlations and physically is 
explained by the well-known sweeping effects, see, e.g., p7| . The RG treatment of this problem 
has been given in Ref. |24] (see also Ref. f8|) and it is readily generalized to our case. It is 
now generally accepted that the intermittency phenomenon leads to a singular L dependence of 
the equal-time correlations, see, e.g. Ref. [p8fl . In Ref. p9|, it has been applied to the simple 



example of the so-called rapid-change model of passive scalar advection |ID| in order to confirm the 
singular dependence of the equal-time correlation functions on L and calculate the corresponding 
anomalous exponents within the e expansion; the results obtained are in agreement with the 



previous results obtained using the so-called zero- mode technique The generalization of 

these results to more realistic models like the stochastic Navier-Stokes requires a considerable 
improvement of the existing technique and remains an open problem. 
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Figure 1: Graphs of composite nonlocal operators F 1 = v\ = v[{djVi — diV^diA^V tA^didjViVj 
and F 2 = v[(Avi) A~ 1 didjV (v j giving a leading correction to the infrared form kinetic energy 
spectra of weakly compressible developed turbulence. 









V// 



Figure 2: The correlation function (v'vvv) with the contribution of the nonlocal operator 
F 2 = vKAv^A' 1 did jViVj. The shadowed rectangle denotes an arbitrary one-particle irreducible 
diagram with four external legs. 
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